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Abstract 

We consider a generalization of the Schur process in which a partition evolves 
from the empty partition into an arbitrary fixed final partition. We obtain a double 
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analysis. If we appropriately scale the length of the row, the limiting correlation 
kernel changes from the extended Airy kernel. 
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1 Introduction 



Recently in nonequilibrium stochastic models which belong to the one- dimensional Kardar- 
Parisi-Zhang(KPZ) universality class pj, remarkable progress has been made on the under- 
standing of the scaling function. The exact form of the function is obtained and its relation 
to the random matrix theory is revealed [21 [3l HI [5], [6] . The common mathematical structure 
among many of these models is that certain "correlation functions" can be expressed as 
a determinant of a correlation kernel. This structure makes it possible to exactly analyze 
the function. The process having the structure is called determinantal process. 

The Schur process, which was introduced in [7j, is the typical process of the determi- 
nantal processes. Let A be a partition, i.e. A = (Ai, A2, ■ ■ ■ ) where Aj G {0, 1, 2, ■ ■ ■ } and 
Ai > A2 > ■ ■ ■ . For a set of partitions {X^^\ A^^-*, ■ ■ ■ , A^^^~^^}, the measure is defined as 

(2N-1 \ 
n SA(^.-)/A(2.)(a('^'))s,(2.+i)/„.„(a(^^'+i))j s„4.v_i)(a(^^)), (1.1) 

where SA/^(a'^*)) is the skew Schur function of shape A//i in the variables 
(For the definition of the Schur function, see (12.101) .) Like the equation above, probability 
measures described by the Schur function appear in various fields in statistical mechanics. 
It has been known that the Schur function is useful for studying the problems of many- 
body nonintersecting random walk and Brownian motion due to its determinant structure 
and the relation to the combinatorics of the Young diagram [8], [9l [TOl [11], [121 [131 • fur- 
thermore the weight (II. ip and its variant appear in analyses of various nonequilibrium 
processes such as the polynuclear growth(PNG) model [14], melting in three-dimensional 
crystal [12], random tiling process [16j, asymmetric exclusion process [TT] and so on. 

In order to visualize this model, we explain many-body random walk interpretation of 
the process introduced in [7]. At t = and t = 4N, the random walker labeled i {i = 
1, 2, ■ ■ ■ ) is on the position 1 — z on one dimensional lattice Z = {■ ■ ■ , —2, —1, 0, 1, 2, ■ ■ ■ }. 
We regard X\^^ — i + 1 as the position of the ith walker, at time n. Here A^-"'' is the ith 
element of the partition A*^"^ = {X^^\x'^\--- ,0,0,---). On the other hand, the Schur 
function S)^(„+l)/^^{n){a^'^^) is interpreted as the transition weight of the walkers from the 

position {X\"^ — i + l}i=i,2,--- at time n to {A^-"'*'"^'' — i + l}i=i,2,--- at time n + 1. Note that 
the variables a*^"^ characterize the transition weight at each time step. 

Due to the property aJ"^ > A^"^ - 1 > ■ ■ ■ > A^^ - i + 1 > aS"\ - i > ■ ■ ■ , the walk- 
ers do not intersect. Note that the transition weight (the Schur function) has the Slater 
determinant like structure (12.101) . Thus, this nonintersecting random walk has the free 
fermionic feature and we find that (II. ip describes the (free fermionic) model of nonin- 
tersecting random walkers which depart from the positions (0, —1, —2, ■ ■ ■ ) at time and 
arrive at the same positions at time AN by way of the positions {A^"'' —i + l}j=i,2, -- at time 
n = 1, 2, ■ ■ ■ , AN — 1, as depicted in Fig. 1. 

This type of nonintersecting random walk is closely related to the PNG model. In [TH 
[18] . the authors introduced the same type of the random walk model called the multilayer 
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PNG model in order to discuss the equal-time multipoint correlations of height fluctuations. 
In this model, the position of the first walker A^"^ at time n, (n = 1, ■ ■ ■ ,4A^) represents 
the height at position n and at some fixed time. Note that in order to find the dynamics 
of the first walker, we need the information about the correlation among the first walker 
and the other walkers. Thus, we focus on the correlation function (see (12.41) ) of the Schur 
process in this paper. We are also interested in seeing how the initial or final configurations 
of the walkers affect the dynamical behavior. In the study of the PNG model with an 
external source [20} [21], it has been known that the property of the height fluctuation 
changes if the external source is larger than a critical value. As we discussed in ^21j, the 
external source in the PNG model may be related to the initial or final configurations of 
the nonintersecting random walk model. 

In order to discuss the problem, we investigate the following process in this paper: 

(2N-1 \ 
n ^A(2.-i) /A(2.) {a^^'^)sxi2.+.) ,xm (a(2^'+^)) j s^(,^-,, /^(4^) (a^^^)). (1.2) 

Here the partition ^^^^^ = {fi^^^\ fi^2^\ ■ ■ ■ , fin^\ 0, 0, ■ ■ ■ ) is arbitrary but fixed. Notice 
that the weig ht ([L2D is reduced to ([H]) if we choose /i^^^) = := (0,0, ■■■)• In the 
interpretation of nonintersecting random walk, the partition /j,^^^^ corresponds to the final 
configuration {^^f^^ —i+l}i=i^2,---- Thus, the weight (II. 2p is a natural extension of the Schur 
process. Fig. 2 illustrates the situation. In this paper, we obtain the double integral formula 
of the correlation kernel of the (dynamical) correlation function (see (12.51) ). Furthermore, 
by the use of the formula, we discuss the asymptotic limit of the correlation function. It 
will be shown that if we appropriately scale the partition /i*-^^-*, the scaling limit changes 
from case (II. ip . 

In [7], Okounkov and Reshetikhin originally derived determinant structure of the cor- 
relation function and the integral representation of the correlation kernel in the Schur 
process (II. ip by the method of using the fermion operators. This method was first de- 
veloped in their previous article [22]. In |Tlj, Johansson considered the height fluctuation 
property of the (multilayer) PNG model which is essentially the same as process (11.11) and 
discussed the derivation of the integral formula by using a property of a Toeplitz matrix. 
Furthermore, the edge scaling limit(see Sec. 2.3) of the correlation kernel is also obtained. 
A simple linear algebraic proof of the integral formula along the approach in [T3] and its 
Pfaffian analog were discussed by Borodin and Rains [23j. In the noncolliding Brownian 
motion, which we can regard as the continuum version of Fig.l, the correlation function and 
its asymptotic form have been discussed in the similar situation to Fig. 1 where Brownian 
particles start and end at one point. It is found that the process is closely related to the 
random matrix theory and orthogonal polynomials [121 l2ll [25] . The connection between 
the dynamics of the first particle and Painleve equation is also revealed [261 [27] . 

In case (11.21) . on the other hand, the asymptotics of the correlation function has not 
been discussed yet since an integral representation of the correlation kernel has not been 
obtained. In this paper, we discuss the generalization of the method in [13] in order to 
obtain the integral formula for the correlation kernel in process (II. 2p although we could 
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consider the generalization of the methods in [71 [23]. In the similar situation of noncolliding 
Brownian motion where the particles start from one initial point and converge to fixed 
final points at the end, it is revealed that the process is related to the random matrix 
model with external source and due to the properties of multiple orthogonal polynomials, 
its asymptotic limit can also be obtained [2H1 [21] . Recently a partial differential equation 
describing the dynamics of the first particle is also obtained in this situation [29]. However, 
process (11. 2p includes the parameters a*^*^ by which we can change the transition weight 
every time step. Therefore, (11.21) is a more general process in the sense that it includes 
both a^*-* describing the temporary inhomogeneity and boundary parameter 

This paper is arranged as follows. In Sec. [21 we discuss the background of this study 
especially the result about the process defined by products of determinant [Mj and give 
our main results: the double integral formula of the correlation kernel (Theorem 12. ip and 
its edge scaling limit (Theorem 12. 2p . The double integral formula of the correlation kernel 
is derived in Sec. [31 In Sec. [H we discuss the edge scaling limit of the correlation kernel 
applying the saddle point method to the integral formula. The Concluding remark is given 
in Sec. [51 



2 Correlation function 

2.1 Determinant representation of correlation function 

The determinantal structure of correlation function has been derived for various stochastic 
processes such as the Schur process (II. ip [7] and random matrices [30]. For later discussion, 
we consider a class of general measures which contains the above measures and describe the 
result obtained in [H]. This is obtained by generalizing the discussion about the random 
matrix theory in |31j. 

For the set {x*^*^}j=i,... ,4Ar_i, where x*^*-* = {xi\x2\--- we consider the weight 

defined by the products of determinants, 

4iV-l 

W^M,4iv({x«}.=l,..,4iV-l) = n {^r,r+l{xf\x^;^'^)^^ , (2.1) 

r=0 

where {(j)r,r+iix,y)}r.=i,... ,m are some functions on and we fix x^^^ and x^^'^^ such that 

xf ^ > ^ > ■ ■ ■ > k = 0, AN. (2.2) 

(In [H], some condition is assumed for the function (f)r,r+i{x,y) such that all objects that 
appear in the discussion converge.) 

By using Lindstrom-Gessel-Viennot's method [321 [33] (see also [31]), we find that under 
condition (12. 2p . the weights assigned for all configurations but those satisfying 

X? > X? > ■ ■ ■ > xg (2.3) 
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vanish. Namely, when we interpret x"p as the position of the particle labeled i at time j 
as in the previous section, the weight does not vanish only in the case where the particles 
do not intersect. 



The correlation function of the measure 
is defined as 



(1) 



.(1) 



{m-i) 



^(4A^-l)^ 



ly ^ ' , , , ly ... 



R{x\ 



(1) ^(1) A^N-l) (iN-1), 



1 



/ 4Af--l oo 

n E 



v 



i=l (i) [i) 



WM,4Ar({x^*^}i=l,...,4Af-l), (2.4) 



.a;M=-oo 



where Z = Y[t=i ^ (>)__ WM,m{{x^^-'}i=i,--- is the normalization constant. 



•^1 '"" ^"^M ~ ^ 



In [m], Johansson showed that it can be represented as the determinant. 



R{x\ 



(1) 



„(i) 



(4Af-l) 



J4A^-l)^ 



... ''y ... 'T* ... y 



l<r,s<4N-l,0<jr<krfi<js<ks 



= det [Kir-x^},s,x^:^))^^ 

The correlation kernel K(ri, xi; r2, X2) is expressed as 

K{ri,Xi;r2,X2) = k{ri,Xi;r2,X2) - 0ri,r2(a;i, ^2). 

Here, 

E^+i=-oo ■ ■ ■ E^_i=-oo <Pr,r+l{x, Xr+l) " " " (ps-lA^ s-1, V) , for r < S, 



0, 



for r > s, 



(2.5) 



(2.6) 



(2.7) 



M 



K{ri,Xi;r2,X2) = ^ (j)r,,AN{xi,x[*'^^){A ^)ij0o,r2(4''^ ^2), 

A, = 0o,4^(xf (2.9) 

Hence, we find that in general class of measures, which is described by products of deter- 
minant, the correlation function is represented as the determinant. 
Noticing the Jacobi-Trudi identity [35], 



(2. 



sx/M) = det {hx,-t,,+j^i{a)) , 



(2.10) 



where hk{a) is the fcth complete symmetric polynomial in variables a = (ai, a2,- ■ ■), 

hk{a) = ^ Oil ■ ■ -Oj^, (2.11) 



il<--<ik 
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we easily find that the weight (11 ■2p has the form of (12.11) . i.e., products of determinants 
under the following identification: 



x(°) = (0,-l,-2,---), (2.12) 

^(4:N) ^ {mi,m2, ■ ■ ■ jTUn, -u, -n - 1, ■ ■ ■), (2.13) 

= /if^)-z + l, (2.14) 

xO) = {X[^\ Xf - 1, - 2, ■ ■ ■ ), for J = 1, ■ ■ ■ , 4iV - 1, (2.15) 

\h (r+i) {r)(a*^'"+^)), for r even, 

.*,„«(.r.4*")= /' , ,/ (2.16) 

)' for r odd. 

Note that in (12.101) . the rank of the determinant can be infinite (infinitely many walkers 
can move); hence, in this case, M in (12.11) and (12. 8p is infinity. Thus, under this identifica- 
tion (l212D - fl216D with M ^ oo, the result ([23])-([2lD is apphcable to case 

2.2 Double integral formula 

One of the purposes of this paper is to obtain a double integral formula of the correlation 
kernel in case (II. 2p . which is useful for the analysis of the scaling limit as N goes to infinity. 
In case of the Schur process (II. ip . the integral formula was first obtained in [7] by using the 
fermion operators. In this paper, however, for considering the case (II. 2p . we generalize the 
approach in [H] of calculating the inverse of the matrix A in (12.81) in case (I2.12p -( l2.16p . In 
the case of the Schur process (11.11) . i.e., the case /x = in (12.141) . the matrix ^ is a Toeplitz 
matrix and its inverse can be estimated by using the Wiener- Hopf factorization [Ti] . 
However, in case (II. 2p . where /i*^^^-* is general, we cannot apply the method because A is 
not a Toeplitz matrix anymore, which is the main difficulty of this problem. 

In this paper, we develop the method of estimating A~^ in the case of the weight (II. 2p 
and give the double integral formula for the correlation kernel. In the following theorem, 
we consider the situation 

a« = (aP,a«,---,a«,0,0,0,---), (2.17) 

(i) 

where < < 1 for i = 1, ■ ■ ■ , AN — 1 and j = 1, ■ ■ ■ ,p. We can easily find that the 
conditions above ensures the convergence of the partition function Z in (12. 4p . 
The theorem is summarized as follows. The proof is given in the next section. 

Theorem 2.1. In the case (12321) - (I2A71) andn = 2ui {i = 1,2) m i^M), K{2ui,xi;2u2,X2) 
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and 02ui,2m2(^1) ^2 ) in the correlation kernel {\2.6\i become 
K{2ui,xi; 2m2,X2) 



U^2 -TT 
11 



2^1 



+ — TT E E J2^-iy^'4-'s,uyMa) , (2.18) 

, ^ ' j=i £'=0 b=i / 

02„.,..(x„a:2) = — y ;TT^n 11 , _ ,(2m)^ 11 TT^^' ^'"^^^ 

where Cr^ denotes the contour with radius satisfying ri > r2, hm{o) is the mth complete 
symmetric polynomial (12. lip in variables 

a = (a[ \ \ a\^ \ ■ ■ ■ , Qi \ \ ci2 \ ci2 \ ■ ■ ■ , \ ■ ■ ■ , ap \ \ \ ■ ■ ■ , •* ) , 

(2.20) 

and 

^ (^(4^), ^(4^), . . . , ^(4A^)) = (^^, m2 + 1, ■ ■ ■ , m„ + n - 1), (2.21) 

= (/^r + 1, ■ ■ ■ , /^f-T + 1, /iJ-tT, ■ ■ ■ , /^r^), (2.22) 
^^'^ = ( l,lr--,l ,p,--,q)- (2.23) 

6-1 n-b+l 

As was discussed in [HI [23] , the condition (I2.17P can be relaxed such that the product 

[[m=iY^ (i_„(2fc)^^)]^"2 "^_^(2^-i)^^^) converges. Note that m case rm = 0(the Schur 

process), the second term in (12.181) vanishes and it reduces to the integral formula in |7j. 



m, — 1 



2.3 Edge scaling limit 

The benefit of the representation (I2.18P is that we can take an asymptotic limit of the 
correlation kernel. In this paper, we consider the edge scaling limit of the correlation 
kernel in the special situation 

x(^^) = (m,-l,-2,..-), (2.24) 
a« = ... = a(^^) = (a,0,0,..-). (2.25) 

In (I2.24P with m = 0, the scaling limit of the correlation function was analyzed in [H] 
in order to study the height fluctuation of the PNG model. We discuss how the scaling 
behavior depends on the value of m. 



7 



It has been known that in the case m = 0, the trace of the first particle xf\ t 
0, • • ■ , AN, behaves asymptotically like 



hm -i- = A{t) := ^ + -^V^t/N-it/Ny. (2.26) 

This represents a semi-circle centered at t = 2N. We focus our attention on the fluctuation 
of the walkers' positions Xi, i = 1,2, ■ ■ ■ , around the above limiting value at t = 2N. The 
scaling limit is called the edge scaling limit. Precisely, Xi and 2ui in (12.181) and (12.191) are 
scaled as follows: 

Xi = A{2ui)N + DN'/^^i, (2.27) 
2ui = 2N + 2CN^/^Ti, (2.28) 

where A{t) is defined in (12.261) and 

^ = T^(l + «)^/^ C=^^±^. (2.29) 

The exponent 1/3 (resp.2/3) in (12.271) (resp. (12.281) ) characterizes the one- dimensional KPZ 
universality [21 El II El E] ■ We also scale m in (OlD as 

m = A{2N)N + BN^'^uj, (2.30) 

where 

2a2/3 

^=(l-a)(l-Fa)V3- (2.31) 

Our result of the scaling limit is summarized as follows. The proof will be given in 
Section IH 

Theorem 2.2. In the situation (12.241) - (12. 3 ip . the correlation kernel has the following 
scaling limit 

lim K{2ui,xi]2u2,X2) — - — 

'/C2(ri,6;r2,6) + /o"f/Ae^(^^+-)Ai(ei - A)Ai(6), n + < 0, 
/C2(ri,6;r2,6) - /o"f/Ae-^(^^+-)Ai(ei + A)Ai(e2) (2.32) 
-Ai(6) exp (-^^ + 6(n + ^)) , n + ^ > 0, 

where D is defined in (I2.29p . P (see (14.41) ) zs a factor which does not contribute the deter- 
minant in (12.51) . anc? 

l-/°^t/i.e--(— )Ai(a + z/)Ai(6 + z/), ri<r2. ^ ^ 
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When m = in fl2.24p . which corresponds to case ([ri])(/i(^^) = 0), only the term 
^2(^1, ^i; T2, ^2) is left since this case corresponds to u ^ —00. The correlation kernel 
A^2(ti,Ci; ■'"2,1^2) is called the extended Airy kernel, which appeared in the random matrix 
theory [361 Ell EHl [39] and the PNG model [ll[T8]. Recently, the correlation kernel (Km 
has also appeared in various fields such as the random matrix with external source (in- 
cluding the noncoUiding Brownian motion with pinned initial or final condition) [211 129] , 
the PNG model with an external source ^01 [20] . statistics [H], and so on. 

3 Integral representation for the correlation kernel 

In this section, we give the proof of Theorem 12.11 First, we represent 4'r,r+i{x,y) fl2.16p 
in terms of the vertex operators defined in (lA.Sp . The analysis of the Schur process using 
the vertex operators were developed in [71 [22]. However, note that in our analysis, we 
adopt the one-particle basis whereas in [3 [22], infinitely-many particle state is chosen as 
a basis. Basic properties of the vertex operators and symmetric functions are summarized 
in Appendix A. 

First, we prove fl2.19p . From fl2.10p and (]A.13p . we have 




for r even, 
for r odd. 



(3.1) 



Here, the parameter s^^^ = {s 




■ ■ ■) is connected to a*--'^ in (11. 2p as 



1 



p 





(3.2) 



By using (13.11) and (]A.9p . we can express 02mi,2m2 (12 -Tp as 



(3.3) 




By applying ( 1A.6P and ( 1A.7P to this equation, we easily get (12.191) . 
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Next, we consider the integral form of K. Eq. (12.81) is written as 
K{2ui, xi, 2u2, X2) 



JCr, ^1 JCr, Z2 ^^^^^ 



(271^)2 Jc., z\ 



X 



00 



(2-)^ ia, .1^- 1, ^2 ^^/^ ^-'^ 

2N-1 2N U2-1 U2 

j=Mi j=ui+l k=0 1=1 

Here, 7(2,5) is defined in (IA.12P and C.^- (z = 1,2) denotes the contour centered at the 
origin anticlockwise with radius satisfying r2 < ri. In the second equality, we use (13. 3p 
and flXejl -dCT). 

Thus, what we have to do is to obtain an explicit form of the inverse of the (semi-infinite) 

matrix A (El]) in case ^TTR - ^M) and calculate ES=i -^i' Aj^^l From (EJl), the 
matrix A in our case (I2.12l) - p.l6l) can be expressed as 



Aij = (x,|A+(s)A_(s')|l-^), (3.5) 



where 



f-l£^£". (3.6) 

I 1 — J, for n + 1 < J, 

2N 2N 

A^is) = l[H^is'^'^-'y), A_(.0 = n^-(^^''^)' (3.7) 

3 = (.«, s^^*), ■ ■ ■ , .(^^-1)), s' = (.(2), ,(4), . . . , ,(4^)). (3.8) 
In order to get A^^, we introduce the following matrix: 

A',^ = {l-j\Az\s')rA-\s)\x.), (3.9) 
where P' is a deformed projection operator 

OD 

r = 5^|i-j)(x,|. (3.10) 

Note that when n = in (13. 6p . P' becomes an ordinary projection operator P defined as 

00 

P = 5^|l-j)(l-j|. (3.11) 

i=i 
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By the use of P' and P, we can express the product A' A as 

iA'A),k = {xk\A+is)A_is')FAz\s')rA-\s)\x,). 
Let us focus on the term 



oo oo 



FAz\s')r = Y.Y1 \l-p){l-p\Az\s')\l-p'){ 
p=i p'=i 



Xp' I . 



(3.12) 



(3.13) 



Noticing 1 — p' < and < j\A_^(s')\i >= for j > i, which immediately follows from 
definition (13.71) . we have 

PA:i(s')r = Az\s')r. (3.14) 



Namely, we can change P into the identity operator in f l3.13p . Thus, we obtain 

iA'A)ik = {xk\A+is)rA-\s)\xi). (3.15) 

Note that when n = 0, one can easily find that P'(= P) can also be changed into the 
identity operator in (13.151) from the similar discussion about (13.141) . Thus, in the case 
n = 0, A' A = 1 or, equivalently. A' is nothing but A~^. 

In the case of general n, we can represent the matrix A' A in terms of Cj_fc(s) and 
dj_k{s) defined in (]A.14p and (]A.15P respectively, i.e., 

c,.u{s) = {j\A4s)\k) = {k\A4s)\j), d,.,{s) = {j\A-\s)\k) = {k\Az\s)\j). (3.16) 

By noting Cj^k = c^j-fe = for j — k < 0, we get 



A'A 



a-2n 



(3.17) 
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where the region after the n + 1th column is equivalent to the unit matrix and 



Ylj=l '^mj-mi('5)Cm,fe+j-l(s), 



1 < i < n, 1 < k < n. 



j=n+l ^rrik+j-l 



[s)di^j{s) + XlJ=i c?mj+i-i(s)cmfc+i-i(s), n + 1 <i, 1 < k < n. 

(3.18) 
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Here, is defined in (12.141) . Note that the matrix A' A is almost the unit matrix. We 
easily find that the inverse of the matrix A' A has the following form: 



/ &11 

hi 



B := (A'A)-^ = f 
bo 



bi2 

&22 



\ 



bnl bn2 ■ ■ ■ bnn 

bn+U bn+l2 ■ ■ ■ bn+ln &0 

V ; ; ; 



bo 



(3.19) 



By introducing 



A" := {aij)^j=i, -^"(ji) '■ ji cofactor of A", 



(3.20) 
(3.21) 



we find that the elements bo and bij can be expressed as 
bo = det A" = det C, 



(3.22) 

1 < i < n, 1 < j < n 



( 



det 



an ai2 



0'n2 



a-in ^ 



- YTk=i (^ikbkj, n+l<i, l<j<n 



(3.23) 



Thus the inverse of the matrix A (12. 9p is represented as 



A 



-1 



BA' 
( 



bo 



fci 



bo 



J_ v^fi h A' 

bo l^k=l '^nk-^kl 
X]fc=l bn+lkA'j^i + A!j^^i^ 
X]fe=l bn+2kA'^i + •^n+2, 



J_ V^" J, /[I 
bo l^k=l "nk-^k2 

bi'^ X]fc=l bn+lkA'^2 + •^n+1,2 

fe;^ X]fc=l ^n+2fc-^fc2 + •^n+2,2 



(3.24) 
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By using this expression for A^^, we can represent J2Ti=i ^i' •^ij'^i i'^ p.4p as 



oo oo _ n oo 



i,j=l i=n+l j=l ^ k=l i=n+l j=l 

n n oo 

+ (3.25) 

° A;=l i=l jr = l 

After quite lengthy calculations of this equation, we obtain the following expression 

«,j=l \fc=l j=l £'=0 6=1 / 

X l[l-\z^\s^-'^-'^)^-\z2,s'^'^'>), (3.26) 
i=i 

where C(bj) is the 6j cofactor of matrix C. The proof of this equation is given in Appendix 
B. By noticing (1A.13P and (]A.14p . we can rewrite Q'(s), detC, and C(bj) as 

ce'is) = he,{a), (3.27) 
det C = s^(4jv) (a), (3.28) 

C(6,) = (-1)^+V)M^)(«)' (3-29) 

where a, and u^^^ are defined in fl2.20p - fl2.23p . respectively. 

By substituting flX^ - flX^ into (IS3D and noting and f[02|) . we get the desired 
expression fl2.18p . 

4 Asymptotic analysis 

In this section, we give the proof of Theorem 12.21 We investigate the asymptotics of the 
correlation kernel by the use of the saddle point analysis. 

In case fl2.24p and fl2.25p . the correlation kernel defined in (12.180 and (12.190 reduces to 

K{2ui,Xi;2u2,X2) = k{2ui, Xi;2u2, X2) - (I)2m,2u2ixi, X2) 

1 f dzi f dz2 zi (l-a/zi)"i(l-az2)^^""' 



{27iiY Jc^^ z\^''^ Jc^^ 4-^^ z, -Z2{1- az,Y^--^{l - a/z2) 



U2 



+ — f il-azY^-^Hl~alzY'^ 
2ni Jc, ^ J K / J 



^Kn-A^J_ f dzi jl ~a/z^r- J_ f dz2 (1 - a;^2)^^-"^ . . 
^ hUa^^) 27rz Jc, zl^^'^' (1 - az^^^-^ 2m Jc, z^-^' (1 - a/z2Y^ ' ^ ' ' 



13 



where Cr denotes the contour enclosing the origin anticlockwise with radius r and hj{a'^^) 
is the jth complete symmetric polynomial (12. lip in variables 



= ... (4.2) 



2N 



In [m], it was shown that the first two terms of the correlation kernel above converge 
to the extended Airy kernel in the edge scaling limits (12.271) and (12.281) . 

DN^/^ f dzi f dz2 zi (1 -a/2i)"i(l -a^a)^^"''" 

lim 



N-.00 (27r2)2 Jc^^ ^l+^i Jc^^ ^1""' ^1 - ^2 (1 - azi)2^-«i(l - a/z2r^ 
DN^/^ r dz 

+ / ^TT^(l-«.)---(l-a/.)— = Px/C,(r„ei;r„6), (4.3) 

where /C2(ri, ^i; r2, ^2) is defined in (12.331) and P is a prefactor which does not contribute 
a determinant, 

2(l+a)2/3jv2/3(^ ) 

P = (1 - a) ^^73 g(rf-r|)/3+C2r2-€in_ (4^4) 

Thus, we only consider the asymptotics of the last term, 

j=l 

where 

1 f dzi (1- alzxY^ , 



^ /• ^^2 (l-a^2)'^~"^ 

under the scaling (I2.27p - (l2.3ip . We scale j as 

J = DN^^^X. (4.8) 

Here, D is defined in (I2.29p . This scaling makes a dominant contribution to the asymptotic 
limit of fl43D . 

First, we consider the function tpi{xi — j) . By substituting (12.281) into (14. 6p . we have 



■ft AT 



exp(iV/^,^(;.i)), (4.9) 



where 



f^Jz) = {1 + P) \og{z -a)-{l-P) log(l -az)-{fi+l + P) logz, 

(3 = CN~^/^Ti, (4.10) 
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and we fix /i as 

/i = + VT^) = ^(2ni), (4.11) 

1 — 

where A{2ui) is defined in ((226]). In (l410ll . C is defined in (E^l). We find that due 
to (14.111) . two saddle points of the function fi3,^{z) merge to the double saddle point Zc{P), 



zM = f±g+"^. (4.12) 



We also have the relations 



dz dz"^ ' dz^ zl{l3) 

Since the main contribution to the integral in (14. 9 p is given around z ~ Zc{l3), z\ may be 
transformed to 

Then, we obtain the relations 

/ / fi'i (z (3)) 

exp(iV/^,^(zi)) ~ exp j^N l^f^A^^ + ^'^ g " 

= expiNfpAzciP))) exp (^^wl^ , (4.15) 
^ exp{{^,-X){tWi-n)). (4.16) 



xi-j+l-^tAf 
^1 



Zl 

In due]), we used ([MZD- By using fimjl - (061) . and 

exp (iV/^,^(z,(/3))) ~ (1 - a) ^T7^ exp(r3/3), (4.17) 

we get 

DiV^/=^Vi(a;i - j) ~ (1 - a) exp - + nX) x Ai(^i - A). (4.18) 
Similarly ip2{x2) becomes 

DN^/^ij2{x2) ~ (1 - a) ^^73 exp (-r|/3 + ^7-2) x ^(^2). (4.19) 
Here, in the equations above, we used the integral representation of the Airy function 

/oo 

^^gixA+|A3_ ^4 20) 
-00 
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Next, we also apply the saddle point analysis to hm-j{ci ) / hm{oL ) in (14.51) under 
the scaling fl230|) and dS]) as ^ oo . From (1X121) . flAlOl) . and (IXSOD . is 
expressed as 

/i .fc^27VN ^ [ -Ng(z) 

where the contour Ci represents the unit circle surrounding the origin anticlockwise, and 
giz) = 21og(l - az) + + \ogz, 5 = (4.22) 

The saddle point Wc of g(z) is 



2a + (l-a)S , 

= 9 ^ n ^- (4-23) 

/a + a(l — a)o 



When we deform the path of z to 

1 



, 1 — a;\ cij 1 — a 

1 + .,1/2 An /2 ^^ ~ 1 + 7T^ + .,i/2/vn/2 ^^> (4-24) 



we find 



exp(-iV(7(2;)) ~exp(-iV(7(^,)-u;2), ~ (4.25) 



Therefore, the asymptotic form of hm-j{oP'^) becomes 



h^_j{a^^) ~ — ^ (l-a)-2^exp ^ + c^VS + cjA . (4.26) 



Thus, we get 

r;^" ^ -exp(Ac.). (4.27) 
From dHHD, dHHD, and KTJ\ . we obtain 

D^'" E -FtW^i^^i - ^>2(a:2) ~ P X / dAe^(-+-)Ai(ei - A)Ai(6). (4.28) 

■ ^^ iT'm[0^ ) Jo 

Here the prefactor P is defined in (14. 4p . From (14. 3 p and (I4.28P and noting that P does not 
contribute to determinant calculation, we finally obtain the desired expression. However, 
it is obvious that this expression is valid only for the case Ti + u < 0{i = 1, 2) because in 
the case Tj + > 0, the integration in (14.280 is divergent. 
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In order to get the asymptotic form in the latter case, we need another expression 
for (14.51) . By using (1A.19|) . we rewrite the term J2^=i hm-j{a'^^)z~^ in (14.51) . 

h^.,{a'^)z^ = _ - Y: Kn-A»^^)z^- (4.29) 

j=l ^ ' ' j=—oo 

Thus, we obtain 

(4.30) 

where 

^(^i) = ^ £ ;t^(1 - o.lz,r--\l - az,r--\ (4.31) 

and '?/'i(x)(resp. '?/'2(a^)) is defined in (14. 6 p (resp. (14.71) ). 

We can easily get the asymptotic forms of v?(x) by using the saddle point analysis in a 
manner similar to the derivation of (I4.26p . The result is 

;i - 2NY^^{\ - a) ^^73 exp — ^ — - - rluj - r^u? + 6^ • 



27rV2 a ^ / V 7 (l + a)2/3 

(4.32) 

From (ICTj) . (giSSD, and (gJl, we have 
Analogous to the derivation of (I4.28p . we also have 

-DN'" E r; 2^^^ 1(3:1 - j>2(x2) ~ -p X / rfAe-^(-^+-)Ai(ei + A)Ai(e2). 

i^myc^ ) Jo 

(4.34) 

From these two relations, we get the desired expression also in the case + > 0. 



5 Conclusion 

In this paper, we have studied the correlation function (12.40 of the process (II. 2p . This 
process is the generalization of the Schur process (II. ip [3| in the sense that in the picture 
of nonintersecting random walk, the walkers end at fixed sites, as depicted in Fig. 2. At 
first, we have obtained the integral representation of the correlation kernel (Theorem 12. ip . 
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The result is the generahzation of the one in the Schur process obtained in [TJ [TU [23] . The 
derivation is given in Sec. [31 The key technique is the calculation of the inverse of the matrix 
A (I2.9P with semi-infinite rank. This technique can be regarded as the generalization of 
the one discussed in [I^. Next, by using the integral representation, we have obtained 
the edge scaling limit of the correlation kernel in the special case fl2.24p and (12.251) . The 
result is summarized as Theorem 12.21 and the proof is given in Sec. [H We have found that 
the limiting correlation kernel is equivalent to the one obtained in other determinantal 
processes such as the random matrix with external source, the PNG model, and statistics. 
We list the future problems as follows. 

1. In this paper, we have discussed the asymptotics of the correlation kernel only in 
the simple case (I2.24p and ( ]2.25p . although the process ( 11. 2p has many parameters 
yi^(4Af) g^j^j {a^*^}j=i,... ,4Ar. It would be interesting to investigate the limiting behavior 
of the correlation kernel in a more general situation and how it depends on these 
parameters. 

2. The correlation kernel may be closely related to the orthogonal polynomials. In the 
Schur process, the correlation kernel (the first term in (I2.18P ) can be represented in 
terms of the Meixner polynomial in the single time case [12]. On the other hand, it 
has been recently revealed that the correlation kernel corresponding to the random 
Hermitian matrix with external source can be expressed in terms of the multiple 
orthogonal polynomials [1^. In the discretized process (II. 2p . is the correlation kernel 
related to any discrete analog of the multiple orthogonal polynomial? 

3. In Theorem 2.1, the second term in (12.180 . which is due to the partition /i*^^^) in (II. 2p . 
can be expressed by using the Schur function. This fact raises our hope that there 
exists some deep connection between the process (II. 2p and the theory of integrable 
system and representation theory. In particular, it is interesting to view the problem 
in perspective of the Kadomtsev-Petviashvili and Toda Lattice hierarchies [1^ [15] . 
Recently, the relationship between the stochastic processes such as the random turn 
walk and soliton theory has been discussed in [45j. This approach may be useful for 
studying this topic. 



A Vertex operators and symmetric functions 

Let Cj and Cj (i, j G Z = (■ ■ ■ , —1, 0, 1, ■ ■ ■ )) be the creation and annihilation operators 
that satisfy the fermion anticommutation relation 



{Cj, cj} := Cic] + c\ci = Sij, (A.l) 
One particle state \i) is defined as 



{q,c,} = {c1,c;} = 0. (A.2) 



\i)=c\\n), (A.3) 
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where \ fl) is the vacuum state without particle. 

The mode operators c{z), c^{z) and the vertex operators H±{s), where s = {si, S2, ■ ■ ■), 
are defined as follows: 



00 00 



cUz 



i=— 00 i=— 00 

(00 \ 00 

YsnP±n\ , where /3±„ = 4±nCfc, (^ = 1, 2, ■■■). (A.5) 
n=l / fc=— 00 

We summarize the basic properties of these operators, which are used for our discussion, 

H±ism = \n), (A.6) 
H±{s)c^{z) = ^{z^\ s)c\z)H^{s), (A.7) 



00 



H±{s)c{z) = j{z^\ s) ^c{z)H±{s), where 7(2;, s) = exp ^ , (A. 



,n=l 



In ( ]A.9[) . the symbol = means an equality on one particle space. Note that (1X91) IS 
different from the ordinary one. As discussed in [7j when the product H^{s)H-{s') appears 
in another vacuum state where infinitely many particles are occupied up to the origin, the 
following relation holds: 

H+{s)H^{s') = e^-i (s')i/+(s). (A.IO) 

When we set the variable s of the vertex operators, 

n 



we find that 7(z, s) in (lA.Sp is given as the generating function of the complete symmetric 
function hj{a), 

n 

-L (Xj Z 



1=1 



By using this and the properties flA.6p - flA.8p . we can describe hj{a) as 

h^.,{a) = {j\H4s)\t) = {t\H_{s)\j), (A.13) 

under the parameterization flA.lip . 

We also often use the following functions: 

c, .,{v) = {j\A+{v)\k) = {k\A4v)\j), (A.14) 

d, ^,{v) = {j\A-^\v)\k) = {k\Az\v)\j), (A.15) 
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where 

2N 



A±(t;) = n^±(^^'^ (A.16) 
(2^)). (A.17) 



The properties that we often use are summarized as follows: 

E'^ / N , / N 1 1, for A; = 0, . . x 

c,{v)d,.,{v) = {' , (A.18) 

oo 2Af oo 2Af 

Y,c,[v)z^ = l[^{z,v^'^), Y.'^M^' = W^-^zM""). (A.19) 

1=0 i=l i=0 2=1 



B Proof of (B2ED 



In this appendix, we derive (13.261) by deforming the right hand side of (13.251) . By noticing 

oo 

fc=i 

and ( 1A.19I) . we rewrite the first term of (13.25^ as 

cx> oo oo 2A'' oo 

i=n+l j=l k=l j=l i=n+l 

oo ^ X fc 2Af 

k=n+l ^^1/ j,=l 

n 00 27V 

+ Y4 E 4.+.-i(^);^rn^~'(^2,^(^^')). (B.2) 

fc=l i=n+l j=l 



We need a quite lengthy calculation for the second term of (13.251) . At first, by using (I3.18P 

00 

i=n+l ■ 



and (I3.23p . the part X^i^n+i ^ term can be represented as 



i=n+l j=l i=n+l 

n 00 /i—n—1 n 

= -E^J'^ E I E Cr„,+i-l-Ks)df(s) +^C„^+a-l(s)d„„+i_l(s) j 2;f 
j=l i=n+l \ £=0 a=l 



(B.3) 
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In this equation, we notice that the term Xli^n+i ^£=0 ^ Cmj+i-i-e{s)de{s)zi * becomes 



E E 

i=n+l e=o 



=n+l 1^1=0 ^-nij+i- 

oo i~n-l rrij+n-l 2N 

c„^,+i_i_£(s)d^(s)2;f* = - ^ c„.+^„_i_£(s)zf-" Jj7-^(2;f \ s^^^"^)) + z^'~^ . 

(B.4) 

The derivation is given as follows. By noticing ( 1A.18I) and (1A.19I) . one gets 

oo i— n— 1 

$Z X] c,„,+i-i-K'S)c?£(-s)2;r 

oo nij+i—l oo mj+n-l 

= - ^ 5Z Cm^+i-l-iis)deis)z];' = - ^ ^ C„^+n-l-£'('5)c^£'+i-n(-5)V 
j=n+l e=i-n i=n+l £'=0 

m,j+n-l / oo 



£'=0 



\i=n+l 
'2N 



I' 



J2 Orn,^n-l-,{s)z('- \{r\z^\ 3^'^-'^) -J2d,{s)z^^ 



l"=0 \i=l 



(B.5) 



i=0 



By noting again ( lA.lSp . we find 

nij+n—l £' 



(B.6) 



£'=0 



i=0 



Then, we eventually obtain ( lB.4p . On the other hand, from ( IB. II) . the other term A'j. 
in the second term in f l3.25p is 



2N 



3=1 



6=1 i=l 



(B.7) 



Thus, from (|R3|) . (iRil) . and flBTD . we get 

^ n oo oo 

E -r^^^E^^.^^ 

A;=l j=7i+l j=l 

^ n n oo 

—j^EE^-^r'-'E^; 

'J k=l j=l c=l 
^ n n oo n A; 2^ 

^EE^j-^ E E^"^^+"-i'^^)^™«+^-i(^)^i~'E^™''-™fe'^^)^2n^~^'^^2,s(^*)) 



k=l j=l i=n+l a=l 
n n / nij+n-l 



fe=l i=l 
k \ 2N 



+ rEE ^^-^ E Cm,+n-l{s)z[--Yd„^{s)z', ] \[^-\z^\ S^'^-'~^)r\z2. 



k=l j=l 



b=l 



1=1 



{B.8] 
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Note that the first term in this equation and the last one in fl3.25p cancel each other. We 
can also find that the second term in (]B.8|) cancels out the second one in (]B.2|) by the 
following discussion. By using the properties 



^^^^ 



■m,,-mk 



[s] = C, 



fc=i 



^Cm^.+a„l(s)C(fcj) = (5afedetC = 5abho, 



(B.9) 
(B.IO) 



where C(bj) is the hj cofactor of the matrix C (13.211) . we have 

^ n n oo ra k 

k=l j=l i=n+l a=l 

n oo 2N 



27V 



i=l 



(B.ll) 



a=l i=n+l i=l 

Furthermore, from relations flB.9p and fIB.lOp . we can also rewrite the third term in fIB.Sp 



as 



mj+n—l 



i=l 



° fc=l j=l V £=0 fe=l 

n n " / \ ' 

i5:5:c.w--"E'?<..4+E(f) 

2N 

^X{l-\z-,\s^'^-'^)r\z,,s^''^). (B.12) 

i=l 

From flRSj) . flRTD) . and flRT2|) . we get 



^ 71 oo OO 



kjZ2 



k=l i=n+l j=l 

n n oo 



n n oo n oo 2A'' 

rEE^.^^r^"E^^e-2-E^2 E rf™.+.-i(^)%-^n^"(^2,.(^^)) 



fc=i i=i 

m, — 1 



c=l 



a=l j=n.+l 



i=l 



X 



j=l £'=0 6=1 6=1 ^ 

n7-^(^r\^(^^-^))7-^(^2,.(^^)). 



(B.13) 



i=l 



Hence, from (jS^SD, and (iRlSl) . we finally obtain ([326]). 
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Figure Captions 



Fig. 1: The random walk interpretation of the weight (11.11) . Note that the ith walker 
starts at the point + 1 and returns to the initial point at the end by way of the 
position A-"^ — i + 1 at time n. In this example, the partitions A'-"^ {n = 1,2,3) are 
A(^) = (4, 4, 2), A(2) = (2, 1, 1), and A^^) = (4, 4, 3, 2, 1), respectively. 

Fig. 2: The random walk interpretation of the weight (11.21) . In this case, the ith walker 
starts at the point 1 — i and ends at the point nf''^^ —i + 1, which is arbitrary but fixed. 
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